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The  light-harvesting  and  charge-transfer  abilities  of 
Photosystem  1  (PS1)  have  generated  interest  in  the  development 
of  this  system  for  alternative  energy  production  and  energy 
conversion.  We  describe  multi-scale  computational  approaches 
that  were  used  to  study  electron  transfer  at  the  PSI-biological 
and  inorganic  interfaces  and  to  model  this  large  protein 
complex.  Our  work  in  the  development  of  an  all  molecular 
dynamics  model  of  the  PSI  monomer  is  shown  and  compared 
to  the  published  experimental  and  ONIOM  optimized  models 
with  differences  noted  in  protein  and  ligand  structure,  electron 
branch  characterization,  and  ionization  and  orbital  potentials 
in  the  P700  chlorophylls.  Differences  between  the  docking 
of  cytochrome  C6  and  plastocyanin  to  PSI  using  established 
docking  algorithms  and  molecular  dynamics  are  described. 
Finally,  dipole  calculations,  luminal  surface  hydropathicity  and 
polarity  characterization  were  used  to  predict  improvements  in 
surface-assembled  monolayer  design. 

Keywords:  Photosystem  I;  cytochrome  C6;  computational 
chemistry;  molecular  modeling;  docking;  quantum  mechanics; 
self-assembled  monolayer;  photoelectrochemistry;  electron 
transfer 


©  2013  American  Chemical  Society 

In  Applications  of  Molecular  Modeling  to  Challenges  in  Clean  Energy;  Fitzgerald,  G.,  et  al.; 
ACS  Symposium  Series;  American  Chemical  Society:  Washington,  DC,  2013. 


Introduction 


Oxygenic  photosynthesis  provides  the  energy  source,  directly  or  indirectly, 
for  all  complex,  multi-cellular  organisms  on  Earth.  This  biological  process 
couples  the  energy  derived  from  photo-excitation  to  split  water  and  reduce  CO2 
in  order  to  create  complex  carbohydrates  in  plants,  algae,  and  cyanobacteria. 
Two  large  transmembrane  protein  complexes,  photosystems  (PS)  1  and  II,  are 
instrumental  in  the  first  steps  involved  in  charge  separation  and  the  resulting 
translocation  of  electrons  across  the  thylakoid  membrane.  The  first  step  of  the 
light-driven  reactions  in  photosynthesis  occurs  in  PSII,  where  light  energy  is 
harnessed  by  the  chlorophylls  within  PSII,  followed  by  the  oxidation  of  water 
to  produce  protons,  oxygen  gas,  and  electrons.  The  electrons  produced  in 
this  reaction  are  quickly  transferred  through  a  chain  of  cofactors  and  carriers 
(including  plastoquinones,  plastocyanin,  and  cytochromes)  to  the  other  light 
harvesting  complex,  PSI.  Depending  upon  the  organism,  PSI  receives  the 
transferred  electrons  from  either  plastocyanin  or  cytochrome  C6  at  the  P700 
chlorophylls  on  its  luminal  side  and  rapidly  transports  the  electrons  to  the  terminal 
iron-sulfur  cluster,  FB,  at  its  stromal  side.  There  the  electrons  are  used  to  reduce 
NADP+  to  NADPH  in  the  ferrodoxin-NADP+  oxidoreductase  complex.  The 
transmembrane  potential  created  from  the  electron  transport  through  the  two 
photosystems  and  the  proton  concentration  gradient  is  further  used  to  drive  ATP 
synthesis  and  other  cellular  processes. 

PSI  has  become  popular  in  the  solar  energy  conversion  literature  recently 
because  of  its  light  harvesting  ability  and  charge  transfer  properties.  Most 
importantly,  PSI  is  able  to  transport  electrons  across  the  thylakoid  membrane 
with  an  internal  quantum  efficiency  near  unity  (7).  Although  the  comparison  of 
energy  conversion  efficiency  among  disparate  processes  is  not  not  straightforward 
(2),  the  harnessing  and  enhancement  of  natural  photosynthetic  processes  is  of 
great  interest,  and  the  integration  of  PSI  into  chemical  or  electrical  systems  for 
the  conversion  of  light  energy  to  chemical  energy  has  been  pursued  for  several 
decades.  These  efforts  were  initially  focused  on  the  optimization  of  FE  production 
in  algae  (3-6)  through  metabolic  (7-9)  and  later  genetic  controls  (8,  10).  More 
recent  work  has  focused  on  decoupling  photosynthesis  by  crosslinking  PSI  with 
hydrogenases  (11-16)  or  metal  complexes  (77,  18)  in  order  to  direct  electron 
flow  to  drive  hydrogen  evolution.  Despite  this  increased  interest  in  PSI,  the 
fundamental  understanding  of  the  dynamic  docking  of  transfer  proteins  to  PSI, 
as  well  interactions  between  PSI  and  non-native  environments  is  lacking.  This 
work  focuses  on  elucidating  the  fundamental  interactions  between  PSI  and  its 
surrounding  non-native  environment  in  order  to  improve  the  integration  of  PSI 
into  energy  generating  devices. 

Computational  analysis  of  proteins  and  their  interaction  with  surrounding 
environments  is  highly  valuable  in  predicting  the  dynamic  behavior  of  these 
complex  biological  systems.  Interest  in  Photosystem  I  in  particular  has  increased 
dramatically  with  the  release  of  the  2.5  A  resolution  cyanobacterial  PSI  structure 
in  2001  (79).  The  extremely  large  size  of  the  protein  complex,  approximately  one 
megadalton  for  the  cyanobacterial  trimer,  makes  computational  characterization 
highly  challenging.  Much  of  the  earlier  work  used  the  published  structure  to 
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develop  models  of  spectral  properties  and  energy  transfer  kinetics  (20-22).  Ivashin 
et  al.  optimized  select  ligands  (a-chlorophyll  a,  phylloqionones)  and  residues  at 
the  B3LYP/6-31G*  level  in  combination  with  ZINDO/S  level  electronic  coupling 
calculations  to  study  differences  between  the  A-  and  B-branches  of  the  electron 
transfer  chain  in  PSI  (23).  Much  attention  has  also  been  paid  to  the  advanced 
treatment  of  electrostatics  in  this  system  (24-26).  The  P700  site  has  been  studied 
with  TDDFT  (24,  27)  and  semiempirical  methods  (28).  Lin  and  O’Malley  used 
a  two  layer  (B3LYP/6-31G(d):UFF)  ONIOM  methodology  to  optimize  models 
of  the  phyllosemiqionone  free  radical  in  the  Aia  and  Am  sites  (29).  Canfield  et 
al.  utilized  the  ONIOM  methodology  to  systematically  partition  and  optimize 
the  entire  protein  in  a  sequential  fashion  (30).  This  optimized  structure  was  then 
used  for  a  variety  of  analyses  (31).  Attempts  to  dock  PSI  are  less  numerous  in  the 
literature.  Myshkin  et  al.  used  a  variety  of  standard  docking  algorithms  to  model 
the  plastocyanin/PSI  complex  (32).  Joliet  et  al.  used  a  specialized  algorithm  to 
dock  the  PsaC  subunit  onto  PsaA/PsaB  (ii). 

The  computational  approaches  in  this  study  are  focused  on  understanding 
the  interactions  between  PSI  and  its  charge  carrier  proteins  or  surrounding  non¬ 
native  environment  in  order  to  develop  improved  coupling  (and  electron  transport) 
between  biological  components  (cytochrome  C6/PSI  and  also  between  PSI  and 
physical  components  of  electrochemical  devices).  Simulations  at  multiple  scales 
were  used  to  address  these  challenges.  Modeling  approaches  included  docking, 
all-atom  molecular  dynamics  (MD),  and  quantum  mechanics  (QM)  using  density 
functional  theory  (DFT).  Development  of  coarse  grained  models  when  system 
size  and  sampling  time  scales  exceed  the  current  limits  of  all-atom  MD  are  also 
underway  and  will  be  reported  elsewhere.  Here  we  detail  initial  findings  for  this 
ongoing  project  which  is  the  first  to  report  PSI  dynamics  and  extend  multiscale 
modeling  of  the  PSI  complex  past  the  QM/MM  level. 


Modeling  the  PSI  Complex 

In  cyanobacteria,  photosystem  I  exists  as  a  clover-shaped  trimer 
embedded  in  the  lipid  bilayer  of  the  thylakoid  membrane.  Each  monomer  is 
composed  of  12  protein  chains,  96  chlorophylls,  22  carotenoids  (P-carotene), 
2  phyloquinones,  3  iron-sulfur  clusters,  4  unique  lipid  molecules,  and  a 
multimer-coordinating  calcium  ion  ( 19).  Nine  of  the  protein  subunits  feature 
transmembrane  domains  while  subunits  PsaC,  PsaD,  and  PsaE  compose 
the  stromal  hump.  The  interface  of  the  two  largest  subunits,  PsaA  and 
PsaB,  creates  a  pseudo  C2  axis  at  the  center  of  the  PSI  monomer  complex 
with  pigment  and  lipid  cofactors  arranged  in  two  branches  along  the  axis 
(79).  P-carotene,  the  lipids  (1,2-dipalmitoyl-phosphatidyl-glycerole  and 
1,2-distearoyl-monogalactosyl-diglyceride),  and  most  of  the  axial  chlorophylls 
make  direct  contact  with  PsaA  and  PsaB  (79).  Protein  chains  PsaJ,  PsaK,  PsaL, 
PsaM  and  PsaX  coordinate  antennae  a-chlorophyll  a  molecules  directly  or 
through  bridging  solvent  contacts  (79).  Stabilizing  hydrophobic  contacts  with  the 
carrotenoids  are  made  by  amino  acid  residues  in  chains  PsaF,  Psal,  PsaJ,  PsaL, 
PsaM,  and  PsaK  (79).  PsaF  has  also  been  proposed  to  coordinate  docking  of 
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plastocyanin  and  cytochrome  C6  to  the  P700  domain  found  at  the  center  of  the  C2 
axis  {34).  In  cyanobacteria,  this  protein  chain  is  significantly  truncated  compared 
to  its  plant  complement  {19)  and  interaction  may  be  limited  due  to  its  relatively 
distal  location.  A  graphical  representation  of  the  PSI  monomer  from  both  the  side 
and  stromal  views  is  shown  in  Figures  la  and  lb. 


Figure  1.  Ribbon  representation  of  cyanobacteria l  PSI  (PDB  ID:  1JB0)  from  side 
(la)  and  stromal  views  (lb).  Coloring  was  chosen  to  differentiate  the  protein 
chains  that  compose  the  complex  and  emphasize  the  C2  axis. 
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Initial  PSI  Molecular  Dynamics  (MD)  Model 


An  all-atom  molecular  model  of  PSI  was  built  using  the  high  resolution 
crystal  structure  (2.5  A  resolution)  of  Photosystem  I  from  cyanobacteria  ( 19) 
(PDB  ID:  1JB0)  as  the  initial  starting  trajectory.  Some  modifications  to  the 
structure  were  necessary  and  included:  splitting  the  PsaA  chain  between  residues 
262  and  266  due  to  missing  residues,  inserting  an  alanine  residue  into  position 
33  of  the  PsaK  chain,  splitting  the  PsaK  chain  between  residues  43  and  55  due 
to  missing  residues,  and  removal  of  native  waters  and  calcium  ion  from  the 
model  structure.  Ligand  parameters  and  atomic  charges  were  determined  using 
the  R.E.D.  vIII  (35)  and  antechamber  (36)  programs.  Parameterization  of  the 
iron-sulfur  (4Fe4S)  clusters  assumed  an  oxidized  state  and  atomic  charges  were 
calculated  using  a  tetra-coordinated  methyl  sulfide  capping  groups  whose  group 
atomic  charges  were  constrained  to  zero,  similar  to  work  performed  by  Torres 
et  al.  (37).  RESP  (restrained  electrostatic  potential)  derived  atomic  charges 
(35)  based  on  QM  calculations  were  determined  consistent  with  the  AMBER 
ff99SBildn/GAFF  (39)  force  fields  using  the  HF-6/31G*  basis  set  for  all  ligands 
with  the  exception  of  the  iron-sulfur  cluster.  Due  to  the  known  limitations  of 
the  Hartree-Fock  basis  sets  to  model  metal  centers,  the  4Fe4S  structure  was 
optimized  and  molecular  electrostatic  potential  (MEP)  calculated  using  the  DFT 
BVP86/cc-pVDZ  basis  set.  The  all-atom  molecular  dynamics  simulation  was 
performed  using  the  AMBER  ff99SBildn  force  field  (40)  with  adjustments  to  the 
phi  and  psi  amino  acids  dihedrals  consistent  with  the  Robert  Best  parameter  set 
(41).  The  protein  complex  was  solvated  using  a  surrounding  octahedron  shell  of 
explicit  TIP3P  waters  which  extended  at  least  12  angstroms  from  the  structure  in 
all  directions.  Explicit  salt  ions  using  the  Cheatham  and  Joung  parameter  set  (42) 
were  used  to  neutralize  the  system.  Periodic  boundary  conditions  were  applied 
using  the  particle  mesh  Ewald  method  approach  with  a  small  (1  fs)  time  step  to 
ensure  stability  in  the  large  complex.  Minimization  was  performed  in  two  stages, 
allowing  for  initial  solvent  minimization  in  the  presence  of  restrained  solute 
molecules  and  later  unrestrained  system  minimization.  Equilibration,  similiarly, 
used  three  stages  with  weak  solute  restraints  during  a  slow  heating  stage  to  the 
final  temperature  of  300K  followed  by  500  ps  of  unrestrained  equilibration  to 
remove  any  hot  spots  in  the  simulation  prior  to  production  MD. 

During  the  first  3  ns  of  all-atom  MD  of  the  PSI  complex,  the  all-atom  RMSD 
consistently  increased  from  the  starting  trajectory  (Figure  2).  Comparing  the  final 
snapshot  at  3  ns  to  published  structure,  significant  loss  of  secondary  structure  can 
be  seen  originating  in  the  PsaC  domain  (Figure  2).  Using  scanning  RMSD  (43) 
to  measure  local  fluctuations  on  the  time  course  of  the  simulation,  the  majority 
of  large  motions  found  in  and  near  PsaC  occurred  proximal  to  the  Fb  iron-sulfur 
cluster.  At  this  region  many  of  the  stabilizing  salt  bridges  and  hydrogen  bonds 
between  PsaC  and  nearby  PsaD  and  PsaE  chains  are  formed.  The  large  atomic 
charges  assigned  to  the  iron-sulfur  clusters  create  competitive  electrostatic 
interactions  at  the  protein  tertiary  domains  which  disrupt  neighboring  salt  bridges 
and  eventually  lead  to  structure  deformation.  After  the  results  from  the  QM 
atomic  charge  fitting  were  re-evaluated,  it  was  determined  that  the  molecular 
electrostatic  potential  for  the  iron-sulfur  clusters  was  further  distributed  along 

181 

In  Applications  of  Molecular  Modeling  to  Challenges  in  Clean  Energy;  Fitzgerald,  G.,  et  al.; 

ACS  Symposium  Series;  American  Chemical  Society:  Washington,  DC,  2013. 


attached  thiol  ligands  and  that  the  charge  constraints  resulted  in  artificially  large 
charges  to  the  iron-sulfur  molecules.  Variance  between  ligand-complex  fit  and 
complex-alone  fit  atomic  charges  confirm  the  need  to  include  4Fe4S  bound 
amino  acids  in  the  charge  assignment  protocol.  While  developing  improved 
4Fe4S  parameters,  a  second  model  of  the  PS1  complex  lacking  the  stromal  hump 
(PsaC,  PsaD,  PsaE  protein  chains  and  the  iron-sulfur  clusters)  was  constructed 
to  complete  simulations  focused  on  dynamic  structural  fluctuations  in  the  lower 
portion  of  the  electron  transfer  chain  and  interactions  of  external  moieties  (e.g. 
transfer  proteins  and  SAMs)  with  the  luminal  surface. 


Figure  2.  Variance  of  the  initial  MD  model  from  the  published  structure  shows 
deformation  in  the  PsaC  chain.  In  the  ribbon  diagram  (top)  the  published 
structure  is  shown  in  red,  while  the  MD  model  is  shown  in  white.  On  the  bottom, 
all-atom  RMSD  of  initial  MD  model  from  PSI  X-ray  structure  is  shown  over  the 

simulation  time  course. 
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PSI  MD  Model  Lacking  a  Stromal  Hump 


The  design  of  a  PSI  model  lacking  a  stromal  hump  followed  the  approach 
described  above.  A  one-dimension  all-atom  RMSD  plot  from  the  published 
structure  through  the  68  ns  time  course  can  be  seen  in  Figure  3. 


Figure  3.  All-atom  RMSD  ofMD  model  (lacking  stromal  hump)  from  PSI  X-ray 

structure. 


Several  differences  between  the  computational  model  and  published  structure 
are  observed  in  the  ligand  interactions,  the  electron  transport  paths,  and  at  the 
solvent  interface.  In  general,  the  MD  model  is  a  more  compressed  system 
than  the  X-ray  structure  showing  constriction  of  solvent-exposed  hydrophobic 
surfaces  in  both  protein  and  ligand  molecules.  Several  of  the  most  visible 
changes  are  found  among  the  antennae  a-chlorophyll  a  molecules.  Of  these,  CL1 
1601  (a  chlorophyll  in  the  M  chain)  shows  the  greatest  movement  among  the 
chlorophylls  at  1 3  A  from  its  initial,  isolated  position  to  close  proximity  with  a 
neighboring  chlorophyll,  CL1  1201  (see  Figure  4).  Similar  movements  are  seen 
with  chlorophylls  CL1  1217  and  1302.  CL1  1217  moves  6  A  from  a  staggered 
position  with  CL1  1209  to  bridge  CL1  1209  and  CL1  1218  (Figure  4).  CL1 

1302  moves  8.3  A  to  improve  contacts  with  neighboring  pigments.  The  energetic 
driving  forces  for  other  chlorophyll  movements  are  less  clear:  for  example,  CL1 

1303  moves  7  A  without  visible  improvement  in  local  interactions  with  antennae 
chlorophyll  neighbors.  The  published  structure  of  the  PSI  complex  lacks  many 
coordinates  for  atoms  associated  with  the  chlorophyll  hydrocarbon  tails  and  many 
of  these  positions  were  built  at  the  time  of  model  generation  based  on  a  fairly 
simplistic  energy  minimizer.  It  is  plausible  that  the  80  ns  time  frame  of  this  MD 
simulation  may  not  be  sufficient  to  allow  for  optimal  packing  or  refolding  of  those 
regions  and  this  may  contribute  to  some  of  the  differences  seen. 
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Figure  4.  Localization  of  chlorophyll  molecucles  in  X-ray  (red)  and  molecular 
dynamics  (white)  structures. 


In  Figure  5,  a  side  view  ribbon  representation  of  both  crystal  and  MD 
models  of  PSI  is  shown.  Differences  in  protein  positions  can  be  observed  with 
the  emphasis  on  the  PsaL  chain  which  is  shown  in  dark  grey.  The  movement 
of  solvent-exposed  helices  from  solution  to  the  transmembrane  region  is  fairly 
consistent  with  structural  motions  in  the  MD  model.  The  crystallization  of 
photosystem  I  was  accomplished  using  P-dodecylmaltoside  as  a  detergent  for 
the  lipoprotein  complex.  Detergents  and  surfactants  can  minimize  hydrophobic 
forces  in  protein  molecules  which  may  account  many  of  the  differences  seen 
at  the  solvent  interface.  Furthermore,  some  uncharacteristically  fast  motions  in 
the  solvent-exposed  peptide  transmembrane  helices  are  seen  in  the  MD  model 
which  suggests  that  future  simulations  which  lack  detergent  (such  as  ours)  should 
include  additional  bilayer  lipid  molecules  beyond  those  defined  in  the  X-ray 
structure  to  maintain  ideal,  native  structure. 
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Figure  5.  Structural  differences  between  the  PSI  X-ray  structure  (top)  and  the 

MD  model  (bottom). 
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Electron  transfer  through  the  core  of  PSI  takes  places  through  one  of  two 
electron  transfer  pathways;  branch  A  or  branch  B  (44).  These  branches  begin 
at  the  luminal  pair  of  chlorophyll  molecules  that  are  designated  the  P700  site 
based  on  the  absorbance  maxima  at  700  nm  (7).  From  here,  the  branches 
diverge  to  separate  chlorophyll  molecules  with  decreasing  redox  potentials. 
Along  the  A  branch,  electron  flow  proceeds  from  the  P700  chlorophylls  to 
A_ia  (an  a-chlorophyll  a  molecule),  Aoa  (another  a-chlorophyll  a  molecule), 
Aia  (a  phyloquinone  molecule),  to  Fx,  Fa,  and  Fb  which  are  4Fe4S  iron-sulfur 
clusters  (23).  The  B  branch  uses  similar  intermediates:  A.ib  (an  a-chlorophyll  a 
molecule),  Aob  (an  a-chlorophyll  a  molecule),  Am  (a  phyloquinone  molecule), 
and  convergence  with  the  A  branch  on  the  use  of  Fx,  Fa,  and  Fb(23).  Electron 
transfer  rates  through  the  respective  branches  are  species  specific  and  are  often 
a  function  of  distances  between  intermediates  in  the  branch.  In  Synechococcus 
elongatus,  a  10  fold  increase  in  electron  transfer  rate  (44-46)  is  seen  in  the  B 
branch  compared  to  the  A  branch.  A  variety  of  factors  have  been  investigated  to 
determine  their  role  in  this  behavior,  including  the  contribution  of  surrounding 
residues  on  the  midpoint  potential  ofphylloquinones  Aia  and  Am  (24,  47),  as  well 
as  a  B  branch  specific  inclusion  of  a  tryptophan  residue  (TRP  B673)  which  may 
act  as  an  electron  acceptor  between  Am  and  the  Fx  iron-sulfur  cluster  (19,  23). 
Some  preliminary  measurements  of  the  chlorophyll  and  phylloquinone  distances 
along  the  A  and  B  branches  were  made  to  compare  against  the  published  structure 
(see  Table  I). 


Table  I.  Measured  Distances  between  Electron  Carriers  in  the 
PSI  X-ray  and  MD  Structures 


Measured  Distance 

X-ray  Structure 

MD  model 

P700  chlorophylls 

6.34  A 

5.80  A 

P700a  to  A-ia 

11.72  A 

12.74  A 

A_ia  to  Aoa 

8.13  A 

8.74  A 

Aoa  to  Aia 

8.78  A 

9.89  A 

P700b  to  A-ib 

11.95  A 

12.62  A 

A.  ib  to  Aob 

8.22  A 

7.60  A 

Aob  to  Aib 

8.78  A 

9.37  A 

P700b  to  Aia 

26.0  A 

26.2  A 

While  the  molecular  coordinates  in  the  interior  of  PSI  complex  are  very 
similar  between  the  MD  model  and  X-ray  structure,  some  changes  are  notable. 
In  the  computational  model,  an  increase  is  seen  in  the  distance  between  the  P700 
chlorophylls  and  the  first  set  of  chlorophyll  intermediates,  as  well  as  the  distance 
between  the  second  pair  of  chlorophyll  and  the  phylloquinone  intermediates.  The 
distances  between  the  magnesium  ions  of  the  P700  chlorophyll  pair  has  decreased 
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by  a  half  angstrom,  although  the  chlorophyll  ring  distance  itself  is  relatively 
unchanged  (Figure  6  and  next  section).  The  distance  between  the  B  branch  P700 
chlorophyll  and  the  A  branch  phylloquinone  has  increased  slightly  compared  to 
the  published  structure  and  the  experimental  value  of  25.4  Angstroms  obtained 
from  spin-echo  modulation  experiments  (48).  In  the  X-ray  structure  the  distances 
between  electron  intermediates  is  smallest  in  the  A  branch  while  in  the  MD 
model,  smaller  distances  are  found  in  B  branch  which  may  contribute  to  the 
differences  seen  in  the  electron  transfer  rates. 


Preliminary  PSI  Quantum  Mechanical  (QM)  Model 

Theoretical  studies  of  a-chlorophyll  a  and  similar  porphyrins  have  been  a 
topic  of  interest  for  decades  due  to  their  critical  involvement  in  photosynthesis. 
Numerous  TDDFT  (time  dependent  density  functional  theory)  studies  of  the  single 
a-chlorophyll  a  (Chi  a)  molecule  exist,  which  were  able  to  reproduce  experimental 
adsorption  energies  with  varying  success  (49,  50). 

The  P700  site  in  photosystem  1  is  a  dimer  made  of  Chi  a  and  Chi  a'  (the 
CIO  epimer  of  Chi  a).  Computational  analysis  of  the  P700  dimer  has  also  been 
performed.  It  is  interesting  to  note  that  a  variety  of  P700  models  have  been 
studied.  The  molecular  orbital  analysis  of  Plato  et  al.  (28),  which  was  performed 
at  the  RF1F-INDO  level,  started  from  the  experimental  X-ray  crystal  structure  and 
expanded  to  a  series  of  models  of  varying  sophistication,  although  no  complete 
geometry  optimization  of  the  P700  site  was  attempted.  This  work  demonstrated 
the  effect  of  the  interactions  of  the  chlorophyll  with  the  protein  environment, 
in  particular  Thr  A743,  His  A860,  and  His  B660,  which  they  postulated  to  be 
implicitly  involved  in  the  development  of  spin  and  charge  asymmetry  within  the 
P700  dimer.  Sun  et  al.  (27)  utilized  TDDFT  on  an  optimized  model  of  P700 
obtained  from  the  crystal  structure  and  found  its  excitation  to  be  “intrinsically 
asymmetric”  without  the  inclusion  of  neighboring  residues.  Saito  and  Ishikita 
(57)  used  QM/MM  to  study  the  cationic  state  and  spin  distribution  of  the  P700 
site,  keeping  the  coordinates  of  the  MM  region  fixed  to  the  crystal  structure. 

Gunner  (52)  noted  the  importance  of  accounting  for  the  normal  thermal 
motions  of  photosynthetic  proteins  and  their  subtle  but  far-reaching  effects  on 
redox  potential  and  other  important  properties.  Dreuw  et  al.  (55)  stressed  the 
importance  of  quantum  mechanical  optimization  from  the  experimental  crystal 
structure  when  studying  excited  states  of  pigment-protein  complexes.  The  most 
complete  structural  optimization  of  PSI  was  generated  by  Canfield  et  al.  (30). 
Flowever,  this  work  utilized  the  ONIOM  model  and  optimized  sections  of  the 
complex  in  a  sequential  fashion.  We  are  unaware  of  any  studies  of  the  P700  site, 
or  indeed  any  portion  of  the  PSI  electron  transfer  chain,  starting  from  a  classically 
thermalized  structure  of  the  PSI  complex. 

The  current  work  attempts  to  address  this.  What  follows  is  a  preliminary 
assessment  of  the  physical  and  electronic  structure  of  a  model  P700  dimer  derived 
from  a  2  ns  averaged  structure  of  the  PSI  MD  model  described  in  the  previous 
section.  Results  are  compared  to  similar  data  derived  from  the  experimental  X-ray 
crystal  structure  (1JB0),  and  to  the  optimized  structure  of  Canfield  et  al.  Quantum 
calculations  are  performed  at  the  B3LYP/6-31G(d,p)  level  of  both  the  neutral  and 
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charged  (+1)  P700  complex.  For  ease  of  computation  in  the  single  point  quantum 
runs  only,  the  phytyl  ester  of  both  Chi  a  and  Chi  a'  has  been  replaced  with  a  methyl 
ester,  as  is  common  in  the  literature  (54).  Differences  between  the  thermalized 
P700  structure  and  the  original  crystal  structure  are  moderate  and  can  be  seen  in 
Figure  6. 


Figure  6.  Overlay  of  thermalized  P700  structure  ( standard  coloring  by  element) 
and  X-ray  crystal  structure  (shown  in  blue). 

The  heavy-atom  RMSD  of  our  structure  with  respect  to  the  X-ray  structure 
is  0.61  angstroms.  There  is  obvious  fluctuation  in  sidechains,  as  expected,  and 
some  variation  in  twisting  of  the  5  membered  moities  of  the  chlorin  ring.  The 
most  prominent  difference  is  a  marked  movement  in  plane  of  the  chlorin  central 
magnesium  atoms.  The  Canfield  structure  (not  shown)  has  a  heavy-atom  RMSD 
of  0.29  angstroms  with  respect  to  the  crystal  structure,  also  due  to  sidechain 
fluctuation  and  some  chlorin  ring  twisting.  The  central  magnesiums  atoms  in  this 
instance  remain  slightly  out  of  plane.  In  all  cases  the  average  interplanar  distance 
between  the  chlorin  rings  remains  close  to  the  published  value  of  3.6  angstroms 
(19). 

Ionization  potentials  are  calculated  from  the  electronic  energy  in  the  neutral 
and  cationic  states.  For  all  three  models,  the  calculated  IP  is  within  the  range 
for  variations  in  the  chlorophyll  dimer  P680  in  photosystem  II  calculated  at  the 
B3LYP/6-31G(d)  level  by  Takahashi  et  al.  (55). 

The  role  of  His  A680  and  B660  has  been  discussed  at  length  within  the 
literature  (28,  56,  57).  Within  the  crystal  structure  1JB0  it  is  apparent  that  these 
residues  are  ligands  of  the  P700  Mg  atoms,  as  shown  by  the  short  HisN-Mg 
distances  given  in  Table  II.  Canfield  et  al.  chose  to  honor  this  bonding  pattern  and 
protonate  at  the  delta  position  (30),  while  in  the  current  model  these  residues  are 
protonated  at  the  epsilon  position.  Loss  of  this  explicit  ligand  may  in  part  account 
for  the  relaxation  of  the  magnesium  ions  back  into  the  plane  of  the  chlorin  rings, 
as  well  as  the  larger  HisN-Mg  distance  in  the  current  model  relative  to  the  crystal 
structure  and  Canfield  structure.  This  point  is  of  great  computational  interest,  as 
previous  studies  have  argued  that  complexation  of  the  chlorin  with  the  central 
magnesium  leads  to  the  destabilization  of  the  chlorophyll  a  HOMO  and  LUMO+1, 
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with  little  effect  on  the  HOMO-1  or  LUMO,  leading  to  a  slight  change  in  the 
HOMO/LUMO  gap  (49,  58).  While  the  frontier  orbitals  of  the  current  model 
do  differ  from  the  1JB0  and  Canfield  models,  the  stated  HOMO/LUMO  trend 
is  not  followed,  and  therefore  the  effect  of  Mg  complexation,  while  undeniably 
important,  has  been  lost  in  the  other  effects  from  thermalizing  the  protein 
complex.  It  must  be  pointed  out,  however,  that  the  calculated  HOMO-LUMO 
gaps  of  the  three  models  are  in  very  reasonable  agreement  both  with  each  other 
and  with  previously  published  values.  HOMO-1  and  LUMO+1  values  are  not 
provided  for  the  current  model,  but  differ  only  slightly  from  the  HOMO  and 
LUMO  values,  as  has  frequently  been  noted  within  the  literature  for  Chi  a  and 
similar  systems  (59). 


Table  II.  Physical  and  Electronic  Structure  Properties  of 
Three  P700  Models  at  B3LYP/6-31G(d,p) 


Current  Work 

1JB0 

Canfield 

IP  (eV) 

6.01 

5.78 

5.78 

Frontier  Orbital  energy 
HOMO  (eV) 

-5.11 

-4.88 

-4.89 

Frontier  Orbital  Energy 
LUMO  (eV) 

-2.91 

-2.57 

-2.7 

HOMO-LUMO  eV 

2.19 

2.31 

2.19 

Mg-Mg  distance 
(angstroms) 

5.8 

6.34 

6.26 

HisN-Mg  distance 
(angstroms) 

3.99/4.49 

2.36/2.26 

2.20/2.22 

Heavy-Atom  RMSD  wrt 
1JB0  (angstroms) 

0.61 

0 

0.29 

A  variety  of  additional  measurable  properties  are  intimately  related  to  the 
geometric  and  electronic  structure  of  the  cyanobacterial  P700  moiety  (in  both  its 
singlet  and  triplet  form),  as  well  as  for  additional  sites  within  the  electron  transfer 
chain.  Methods  used  have  included  FTIR  (60-62),  as  well  as  EPR  and  ENDOR 
techniques  (28,  63,  64),  and  some  quantum  chemical  analysis  has  been  performed. 
Work  evaluating  our  thermalized  model  within  the  context  of  this  data  is  underway. 

Biological  Interactions  with  PSI 

In  cyanobacteria,  electron  transfer  to  PSI  occurs  through  the  intermediate 
of  charge  carrier  proteins,  plastocyanin  or  cytochrome  C6.  In  order  to  efficiently 
couple  the  charge  carrier  protein,  cytochrome  C6  in  this  case,  to  PSI  outside  of 
its  natural  environment,  a  detailed  fundamental  understanding  of  this  dynamic 
protein-protein  docking  event  is  vital.  Significant  prior  research  focused  on  the 
binding  of  plastocyanin  and  cytochrome  C6  to  PSI  has  been  published  and  while 
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the  conformation  and  orientation  of  the  plastocyanin-PSI  complex  has  been  well 
elucidated,  the  cytochrome  C6  docking  structure  is  still  unknown.  Ford  et  al.  were 
able  to  characterize  the  binding  site  of  plastocyanin  to  PSI  in  vascular  plants  using 
Fourier  difference  analysis  of  arrays  of  PSI  located  on  the  chloroplast  granula 
((55).  They  propose  that  plastocyanin  binds  close  to  the  center  of  PSI  complex 
between  the  PsaA  and  PsaB  domains  with  a  slight  bias  towards  the  PsaL  subunit 
of  the  complex  (65).  In  this  binding  pocket,  two  solvent-exposed,  hydrophobic 
tryptophan  residues  may  form  an  important  feature  of  the  ligand  recognition  site 
(66).  While  no  structures  of  cytochrome  C6  bound  to  PSI  have  been  resolved  or 
published  in  the  scientific  literature,  cytochrome  C6  has  been  proposed  to  bind 
to  PSI  in  a  similar  manner  as  plastocyanin.  Plastocyanin  and  cytochrome  C6 
share  many  attributes  which  suggest  a  common  docking  orientation  including 
similar  hydrophobic  faces  that  match  the  docking  surface  of  PSI  (67),  similar 
response  to  docking  site  mutations  (68),  and  similar  increases  in  affinity  when 
positive  mutations  are  introduced  to  increase  the  electrostatic  attraction  between 
the  proteins  and  PSI  (68). 

Initial  efforts  to  determine  the  cytochrome  C6-PSI  binding  site  followed  the 
efforts  of  Myshkin  et  al.  (32)  by  using  established  docking  algorithms  to  compare 
the  plastocyanin  and  cytochrome  C6  predicted  docking  complexes.  Molecular 
structures  used  in  docking  (and  later  all-MD  simulation)  were  taken  from  RCSB 
Protein  Data  Bank  and  include  the  crystal  structure  of  Photosystem  1  from 
Synechococcus  elongatus  (PDB  ID:  1JB0)  (79),  an  NMR  structure  of  reduced 
plastocyanin  from  Synechocystic  sp.  PCC  6803  (PDB  ID:  1JXD)  (69),  and  the 
NMR  solution  structure  of  cytochrome  C6  from  Synechococcus  elongatus  (PDB 
ID:  1C6S)  (70).  Due  to  the  limitations  of  the  docking  software,  only  protein 
components  were  used.  Docked  complex  structures  from  the  HADDOCK  (71), 
CLUSPRO  (72),  GRAMM-X  (73),  and  RosettaDOCK  (74)  servers  emphasized  a 
noncannonical,  unlikely  transmembrane  localization  of  the  electron  carriers  with 
respect  to  PSI.  Refinement  of  both  plastocyanin  and  cytochrome  C6  docking  to 
the  luminal  surface  of  PSI  followed  using  a  Python  based  interface  to  the  Rosetta 
docking  program  (75).  These  docking  simulations  implemented  a  3  angstrom 
translation  and  8  degree  rotation  angle  of  the  ligand  molecules,  followed  by  a 
subsequent  alignment  of  the  molecules  to  explore  the  PSI  surface.  4000  docking 
configurations  were  tested  and  the  docking  energy  predicted  by  the  Rosetta 
algorithm  was  calculated  for  each  structure.  This  Rosetta  energy  term  was 
plotted  against  the  calculated  distance  of  the  ligand  ion  (copper  in  plastocyanin 
and  a  heme  iron  in  cytochrome  C6)  and  the  magnesium  ions  found  in  the  P700 
chlorophylls.  Five  representative  structures  were  chosen  that  minimized  both 
terms  and  the  structures  were  compared  as  seen  in  Figures  7  and  8.  In  the  models 
of  plastocyanin  docked  to  PSI,  a  common  docking  face  and  domain  were  seen 
in  all  five  models  which  compared  favorably  with  the  predictions  of  Ford  et  al. 
((55).  In  the  models  using  cytochrome  C6,  at  least  two  different  docking  motifs 
were  seen  which  lacked  a  common  docking  face  and  orientation.  Possible  causes 
for  this  alignment  failure  may  have  resulted  from  several  sources  including:  the 
exclusion  of  metal  ions  which  may  contribute  to  aligning  electrostatic  interactions 
during  the  protein-protein  docking;  the  rigid  docking  preformed  may  not  be  able 
to  capture  structural  changes  that  occur  at  the  time  of  docking  (76)  to  stabilize 
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the  bound  complex,  in  which  case  flexible  docking  should  be  pursued;  or  driving 
forces  in  the  diffusion  controlled  (77)  docking  of  cytochrome  C6  with  PS1  may  be 
sufficiently  small  that  Rosetta  was  unable  to  correctly  estimate  their  value. 


Figure  7.  Combined  structures  showing  docking  convergence  in  the  jive  most 
favorable  docked  PSI-plastocyanin  Rosetta  models. 


To  address  the  possible  causes  of  failure  in  the  rigid  docking  calculation 
of  cytochrome  C6  to  PSI,  follow-up  all-atom  MD  simulations  of  the  docked 
complexes  using  the  most  favorable  five  Rosetta-docked  structures  as  the  initial 
starting  trajectories  were  run.  These  MD  simulations  of  the  docked  structures 
were  completed  using  the  AMBER  modeling  suite  and  Best  phi  and  psi  adjustment 
of  the  ff99S  B  i  I  dn  force  field  as  described  earlier.  Heme  parameters  were  taken 
from  the  contributions  of  Bayly  and  Case  which  are  distributed  with  the  AMBER 
modeling  suite  (36).  Copper  bond,  angle,  and  dihedral  parameters  were  taken 
from  Ross  Walker’s  tutorial  using  poplar  plastocyanin.  After  80  ns  of  production 
MD,  large  differences  were  seen  between  the  plastocyanin  and  cytochrome  C6 
docking  trajectories.  MD  simulations  of  the  docked  plastocyanin-PSI  complex 
showed  little  structural  fluctuations  and  varied  minimally  from  the  starting 
coordinates.  Cytochrome  C6,  on  the  other  hand,  was  seen  to  traverse  nearly  the 
entire  luminal  surface  of  the  PSI  complex.  While  all  the  plastocyanin  docking 
trajectories  maintained  the  original  PSI  contacts,  rotations  in  cytochrome  C6 
molecules  were  evident  with  three  of  the  five  docking  trajectories  adopting  a 
common  docking  orientation  following  the  MD  time  course.  Rotation  of  the 
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cytochrome  C6  molecules  to  the  common  orientation  was  seen  in  two  (of  the 
three)  instances  to  occur  in  proximity  to  the  PsaF  domain  which  may  play  a  role 
in  orienting  cytochrome  C6  prior  to  docking  and  electron  transfer.  Further  studies 
will  be  necessary  to  confirm  this  observation. 


Figure  8.  Rosetta  docked  models  showing  orientation  and  faces  of  the  most 
favorable  PSI-cytochrome  cr,  structures.  Clustering  of  the  docked  structures 
suggest  at  least  two  different  motifs  among  the  predicted  complexes. 


PSI  Surface  Interactions 

Computational  Approaches 

Optimization  of  the  interface  between  PSI  and  inorganic  surfaces  is  necessary 
for  the  effective  electron  transfer  between  the  surface  and  biomolecule.  For 
instance,  attachment  of  PSI  to  gold  nanoparticles  and  unmodified  gold  electrodes 
has  produced  photosensitive  and  photoreactive  systems  (78,  79);  however  this 
direct  attachment  to  the  gold  surface  can  lead  to  protein  denaturation  and  block 
access  to  active  sites.  The  use  of  self-assembled  monolayers  (SAMs)  to  enhance 
binding  and  stability  of  PSI,  as  well  as  control  its  orientation  on  the  electrode 
has  been  shown  in  the  literature  (80,  81),  however  with  limited  fundamental 
understanding  of  the  PSI/SAM  interaction.  Therefore,  a  comprehensive  study  of 
the  surface  assembly  of  PSI  on  various  SAM  modified  electrodes  was  proposed 
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utilizing  experimental  and  computational  approaches  to  characterize  the  surface 
of  PS1  and  to  optimize  SAM  candidates.  Computational  surface  characterization 
was  focused  on  polarity,  luminal  surface  electrostatic  potential,  and  luminal 
surface  hydropathicity. 

Our  previous  analysis  of  the  PSI  complex  structure  showed  a  very 
strong  dipole  oriented  at  slight  tilt  from  the  complex  central  axis,  which  may 
possibly  contribute  to  orientation  of  the  complex  stromal-  or  luminal-side  down 
during  deposition  on  a  surface  (82).  The  Kyte-Doolittle  representation  of  the 
hydropathicity  of  the  photosystem  I  monomer  (1JB0)  shows  distinct  segregation 
of  hydrophobic  residues  (red)  to  the  interior  and  external,  bilayer-exposed  sides 
of  the  complex  (Figure  9).  Flydrophilic  residues  (blue)  congregate  to  both  stromal 
and  luminal  faces  consistent  with  aqueous  exposure.  Notably,  distribution  of  the 
luminal  hydrophobic  residues  (in  red)  is  uniform  across  the  plane. 


Figure  9.  Hydropathicity  analysis  (Kyte-Doolittle  algorithm)  of  side  (10a)  and 
luminal  surfaces  (10b)  of  photosystem  1. 

When  the  surfaces  are  remapped  according  to  residue  type,  most  of  the 
hydrophilic  luminal  regions  are  shown  to  be  composed  of  polar,  non-charged 
residues  (see  Figure  10).  The  P700  docking  site  can  be  found  near  the  center  of 
luminal  face,  composed  of  a  hydrophobic  pocket  surrounded  by  distinct  positive 
and  negative  regions  which  is  consistent  with  predicted  long-range  electrostatic 
steering  of  plastocyanin  to  the  P700  binding  site  (55).  Upon  complex  formation, 
short  range  hydrophobic  and  dipole-dipole  interactions  guide  and  ensure  the 
optimal  configuration  of  the  mediator  prior  to  electron  transfer  (83,  84). 

From  these  surface  mapping  studies,  some  characteristics  for  a  luminal  bound 
probe  or  SAM  were  made  and  later  tested.  The  luminal  surface  of  PSI  is  highly 
polar  suggesting  the  favorability  of  a  charged  or  polar  probe.  The  calculated  dipole 
traverses  the  stromal  hump  showing  a  vector  moving  away  from  the  luminal  face 
which  indicates  that  a  negative  charge  would  be  preferred  over  a  positive  SAM. 
The  fairly  homogenous  distribution  of  hydrophobic  and  hydrophilic  regions  along 
the  luminal  surface  may  indicate  the  benefits  of  an  amphiphilic  molecule. 
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Figure  10.  Luminal  surface  mapping  of  the  PS1  published  structure  (1JB0)  by 
residue  type.  Acidic  residues  are  shown  in  red,  basic  in  blue,  polar  in  green,  and 
nonpolar  in  white.  The  P700  docking  site  is  indicated  by  a  yellow  ellipsoid. 


Experimental  Results 

In  order  to  understand  the  dynamic  effect  of  PS1  surface  functionalities  in 
non-native  environments,  experimental  surface  assembly  studies  were  performed 
on  modified  surfaces.  Alkanethiol  self  assembled  monolayers  (SAMs)  on  gold 
substrates  were  chosen  for  the  surface  modifications  in  this  work  because  of  their 
facile  fabrication,  high  reproducibility  and  easy  of  tailoring  exposed  functional 
groups.  SAMs  were  composed  of  six-carbon  length  alkanethiol  chains  with  varied 
terminal  functional  groups;  amine  (aminohexanethiol,  AHT),  methyl  (hexanethiol, 
HT),  alcohol  (mercaptohexanol,  MHO)  and  carboxylic  acid  (mercaptohexanoic 
acid,  MHA). 
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Gold  substrates  were  formed  via  E-beam  evaporation  of  gold  on  Silicon 
wafers  (3  000 A  with  300A  Ti  binding  layer),  followed  by  gold  island  flattening  by 
intermittent  flaming.  SAMs  were  formed  on  Au  electrodes  by  three  day  assembly 
in  ImM  thiol  ethanolic  solution.  The  SAMs  were  rinsed  with  ethanol,  and  dried 
with  nitrogen  gas.  PS1  assembly  on  the  SAMs  was  conducted  via  90  minute 
adsorption  at  4°C  in  the  dark.  After  adsorption,  all  samples  were  rinsed  with 
deionized  water,  and  dried  with  nitrogen  gas. 

Experimental  analysis  of  PSI  surface  assembly  via  Atomic  Force  Microscopy 
(AFM)  showed  that  PSI  assembled  in  higher  density  to  hydrophilic  surfaces. 
The  hydrophobic  SAM,  hexanethiol  (HT),  showed  negligible  PSI  assembly, 
with  some  PSI  adherence  at  the  gold  island  edge  sites  (Figure  11a).  In  contrast, 
aminohexanethiol  (AHT)  showed  high  PSI  binding  with  approximately  621 
PSI/pm2.  Similarly,  mercaptohexanol  (MHO)  and  mercaptohexanoic  acid 
(MHA),  negatively  charged  hydrophilic  surfaces,  showed  similar  densities  of 
PSI;  626  and  502  PSI/pm2  respectively. 


Figure  11.  AFM  images  of  PSI  adsorption  on  (a)  hexanethiol,  (b) 
aminohexanethiol,  (c)  mercaptohexanoic  acid,  and  (d)  mercaptohexanol  SAMs. 


In  agreement  with  the  computational  based  predictions,  hydrophilic  SAMs 
were  more  effective  than  their  hydrophobic  complement,  however  the  predictions 
of  charge  preference  are  more  complex.  Hydroxyl  (MHO)  and  amine  terminated 
(AHT)  SAMs  show  similar  PSI  adherence  rates  with  a  slight  improvement  for  the 
electronegative  reactive  groups  (in  agreement  with  predictions). 

Mercaptohexananoic  acid  which  is  more  electronegative  than  MHO  showed 
a  decrease  in  PSI  binding  suggesting  that  charged  (or  highly  polar)  probes  may 
experience  electrostatic  repulsion  from  the  PSI  surface.  Due  to  the  limitations  of 
available  SAM  head  groups  we  were  unable  to  test  the  effects  of  amphilicity  on 
PSI  adherence  but  intend  to  address  those  questions  in  future  studies. 
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Conclusions 


Computational  approaches  are  used  to  model  the  photosystem  I  complex 
and  its  interactions  with  biological  electron  carriers  which  can  be  used  predict 
improvement  with  electrochemistry  components.  Long  term  all-atom  simulations 
of  PSI  reveal  contraction  of  solvent  exposed  protein  domains  and  repacking  of 
a-chorophyll  antennae  molecules.  These  changes  were  especially  prevalent  with 
the  PsaL  protein  chain  and  chlorophylls  CL1  1601,  CL1  1217,  CL1  1302,  and 
CL1  1303.  Many  of  these  structural  changes  can  be  attributed  to  differences  in 
simulation  and  crystallization  solutions  and  suggest  the  need  of  lipid  molecules 
or  bilayer  embedding  for  improved  long-term  simulation.  Additionally,  loss  of 
secondary  structure  in  our  initial  PSI  molecular  dynamics  model,  principally 
within  the  stromal  hump,  expose  limitations  in  current  iron-sulfur  cluster 
parameterization  approaches  and  demonstrate  the  need  to  include  full  amino  acid 
ligands  in  the  charge  calculation.  Structure  and  distance  comparison  of  electron 
transfer  branches  show  similarity  to  published  structure  with  a  decrease  in  B 
branch  electron  transfer  distances.  Early  quantum  analysis  of  ionization  potential 
and  orbital  energies  of  the  P700  chlorophylls  from  this  model,  the  published 
X-ray  structure  and  the  Canfield  ONIOM-optimized  model  show  a  consistent 
ionization  potential  and  orbital  energies  with  reasonable  agreement. 

Docking  studies  performed  with  established  algorithms  and  subsequent 
all-atom  MD  simulations  suggest  a  difference  between  plastocyanin  and 
cytochrome  c&  docking  and  interaction  with  the  PSI  luminal  surface.  In  particular 
cytochrome  C6  was  found  to  traverse  and  sample  much  of  the  PSI  luminal  surface 
while  plastocyanin  forms  a  more  static  interaction  at  the  P700  site.  In  addition, 
PsaF  may  play  a  role  in  orientation  correction  of  cytochrome  C6  at  the  luminal 
surface  in  cyanobacteria. 

Computational  PSI  luminal  surface  characterization  and  dipole  calculations 
were  used  to  improve  design  of  surface-assembly  molecules  (SAM).  It  was  shown 
that  SAM  probes  that  were  hydrophilic  and  electronegative  demonstrated  best  PSI 
attachment  to  gold-plated  electrodes;  however,  experimental  verification  indicated 
that  highly  electronegative  or  ionic  head  groups  may  actually  decrease  PSI  counts. 
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